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(FORMULAS OF BRION, LAWRENCE, AND VARCHENKO ON 
RATIONAL GENERATING FUNCTIONS FOR CONES) 

MATTHIAS BECK, CHRISTIAN HAASE, AND FRANK SOTTILE 



Our aim is to illustrate two gems of discrete geometry, namely formulas of Michel Brion [7] 
and of James Lawrence [15] and Alexander N. Varchenko [16], which at first sight seem hard 
to believe, and which — even after some years of studying them — still provoke a slight feeling 
of mystery in us. Let us start with some examples. 

Suppose we would like to list all positive integers. Although there are many, we may list 
them compactly in the form of a generating function: 

(i) x 1 + x 2 + x 3 h — = y^x k 



^— ' 1 — X 

Let us list, in a similar way, all integers less than or equal to 5: 

(2) haT 1 + x° + X 1 +x 2 + x 3 + x A + x 5 = x k = — - — - 

1 - x~ l 

k<5 

Adding the two rational function right-hand sides leads to a miraculous cancellation 



x 5 



tXj Jb iXj Jb Jb Jb 

x 1 — x~ l 1 — xx— 1 1 — X 



(3) — h - = — h — — r = — = X + X 2 + X 3 + X 4 + X 5 . 



This sum of rational functions representing two infinite series collapses into a polynomial 
representing a finite series. This is a one- dimensional instance of a theorem due to Michel 
Brion. We can think of as a function listing the integer points in the ray [l,oo) and 
of (121 as a function listing the integer points in the ray (—00, 5]. The respective rational 
generating functions add up to the polynomial that lists the integer points in the interval 
[1, 5]. Here is a picture of this arithmetic. 
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Let us move up one dimension. Consider the quadrilateral Q with vertices (0,0), (2,0), 
(4,2), and (0,2). 

(0,2), _ # (4,2) 



Q 



(0,0* * *(2,0) 

The analog of the generating functions ([1]) and (T5]) are the generating functions of the cones 
at each vertex generated by the edges at that vertex. For example, the two edges touching 
the origin generate the nonnegative quadrant, which has the generating function 

E *"»* = E*" • E»" = (rbroV 

m,n>0 m>0 n>0 v 1 K y > 

The two edges incident to (0, 2) generate the cone (0, 2) + JR> (0, -2) + M> (4, 0), with the 
generating function 

x m y n = " 

(1 -x)(l -y 1 ) ' 

m>0,n<2 v JK y 1 

The third such vertex cone, at (4,2), is (4,2) +R> (-4, 0) + R> (-2, -2), which has the 
generating function 

4 2 

x y 



(1 - x- 1 )^ - x^y- 1 ) ' 
Finally, the fourth vertex cone is (2, 0) +M> (2, 2) +M> (— 2, 0), with the generating function 

r.2 



x- 



(1 — xy)(l — x^ 1 ) 

Inspired by our one-dimensional example above, we add those four rational functions: 

1 y 2 x 4 y 2 x 2 



(l-x)(l-y) (l-x)(l-y- 1 ) (1 -x- x )(l -x-V 1 ) (1 - xy)(l - x- 1 ) 
= y 2 + xy 2 + x 2 y 2 + x 2 y 2 + x 4 y 2 
+ y + xy + x 2 y + x 3 y 
+ 1 + x +x 2 . 

The sum of rational functions again collapses to a polynomial, which encodes precisely those 
integer points that are contained in the quadrilateral Q. 

Brion's Theorem says that this magic happens for any polytope V in any dimension 
d, provided that V has rational vertices. (More precisely, the edges of V have rational 
directions.) The vertex cone /C v at vertex v is the cone with apex v and generators the edge 
directions emanating from v. The generating function 



x m 



3 



for such a cone is a rational function (again, provided that V has rational vertices). Here 



we abbreviate x m for x\ 



Til ™"12 



Brion's Formula says that the rational functions rep- 



resenting the integer points in each vertex cone sum up to the polynomial 0>(x) encoding 
the integer points in V: 



v a vertex of V 



A second theorem, which shows a similar collapse of generating functions of cones, is due 
(independently) to James Lawrence and to Alexander Varchenko. We illustrate it with the 
example of the quadrilateral Q. Choose a direction vector £ that is not perpendicular to any 
edge of Q, for example we could take £ = (2, 1). Now at each vertex v of Q, we form a (not 
necessarily closed) cone generated by the edge directions m as follows. If w • £ > 0, then we 
take its nonnegative span, and if w ■ £ < 0, we take its negative span. 




For example, the edge directions at the origin are along the positive axes and so this cone 
is again the nonnegative quadrant. At the vertex (2, 0) the edge directions are (—2, 0) and 
(2, 2). The first has negative dot product with £ and the second has positive dot product, and 
so we obtain the half-open cone (2, 0) + M <0 (-2, 0) +M> (2, 2) = (2, 0) +M> (2, 0) +M> (2, 2). 
At the vertex (4, 2) both edge directions have negative dot product with £ and we get the 
open cone (4,2) + M >0 (0,4) + R >0 (2, 2), and at the vertex (0,2) we get the half-open cone 
(0, 2) + M>o(2, 0) + M>o(0, 2). The respective generating functions are 



x 



6„,3 



x"y 



and 



If 



(1 - x)(l - y) ' (1 - x)(l - xy) ' (1 - xy)(l - y) ' " (1 - x)(l - y) 

Now we add them with signs according to the parity of the number of negative (w • £ < 0) 
edge directions w at the vertex. In our example, we obtain 



X" 



x 6 y 3 



y 



x)(l-y) 



x)(l — xy) 



x)(l-y) 



- V 
+ y 
+ 1 



9 9 9 2 2 4 2 

xy + x y + x y +xy 
xy + x 2 y + x 3 y 
x + x 2 . 



This sum of rational functions again collapses to the polynomial that encodes the integer 
points in Q. This should be clear here, for the integer points in the nonnegative quadrant 
are counted with a sign ±, depending upon the cone in which they lie, and these coefficients 
cancel except for the integer points in the polytope Q. 
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The identity illustrated by this example works for any simple polyope — a c/-polytope where 
every vertex meets exactly d edges. Given a simple polytope, choose a direction vector (6R d 
that is not perpendicular to any edge direction. Let E+(£) be the edge directions w at a 
vertex v with w • £ > and E~(£) be those with w • £ < 0. Define the cone 

1C(,v := v + ^ K> w + R <o w • 

This is the analogue of the cones in our previous example. The Lawrence- Varchenko For- 
mula says that adding the rational functions of these cones with appropriate signs gives the 
polynomial <J-p(x) encoding the integer points in V: 

a v (x) = (- 1 ) l ^ (0l ^,v(^) • 

v a vertex of V 

Here, a/c lv {x) is the generating function encoding the integer points in the cone /Q )V . An 
interesting feature of this identity, which also distinguishes it from Brion's Formula, is that 
the power series generating functions have a common region of convergence. Also, it holds 
without any restriction that the polytope be rational. In the general case, the generating 
functions of the cones are holomorphic functions, which we can add, as they have a common 
domain (the common region of convergence). 



Proofs 

Brion's original proof of his formula [7] used the Lefschetz-Riemann-Roch theorem in 
equivariant i^-theory [5] applied to a singular toric variety. Fortunately for us, the remarkable 
formulas of Brion and of Lawrence- Varchenko now have easy proofs, based on counting. 

Let us first consider an example based on the cone JC = M>o(0, 1) + M>o(2, 1). The open 
circles in the picture on the left in Figure [1] represent the semigroup N(0, 1) + N(2, 1), which 




Figure 1 . Tiling a simple cone by translates of its fundamental parallelepiped. 

is a proper subsemigroup of the integer points /C H Z 2 in /C. The picture on the right shows 
how translates of the fundamental half-open parallelepipied V by this subsemigroup cover 
/C. This gives the formula 

a K (x) = a v {x) ■ Y, = (1 _ x , } , 

m,n>0 v /v t " 

as the fundamental parallelepiped V contains two integer points, the origin and the point 
(1,1). 



+ A + /\ = 

A simple rational cone in IR d has the form 

{d \ d 

v + ^ I Ai G R> > = v + Y K> w, , 
i=l J i=l 

where w 1; . . . , G Z d are linearly independent. This cone is tiled by the (Nwi + • • • + Nw^)- 
translates of the half-open parallelepiped 

V := |v + ^A iWi | < A, < 1 

The generating function for V is the polynomial 

a v {x) = x m , 

mevnz d 

and so the generating function for K, is 

a v (x) 



(1 - x Wl ) • • -(1 - x w <0 ' 

aGNwiH hNw d v y v y 

which is a rational function. This formula and its proof do not require that the apex v be 
rational, but only that the generators Wj of the cone be linearly independent vectors in Z d . 
A rational cone K with apex v and generators wi, . . . , w n G Z d has the form 

JC = v + K> wi H h M> w n . 

If there is a vector £ G K d with £ ■ Wj > for i = 1, . . . ,n, then K, is strictly convex. A 
fundamental result on convexity [21 Lemma VIII. 2. 3] is that K, may be decomposed into 
simple cones /Ci, . . . ,/Q having pairwise disjoint interiors, each with apex v and generated 
by d of the generators Wi, . . . , w n of /C. We would like to add the generating functions for 
each cone /Q to obtain the generating function for /C. However, some of the cones may have 
lattice points in common, and some device is needed to treat the subsequent overcounting. 

An elegant way to do this is to avoid the overcounting altogether by translating all the 
cones [5]. We explain this. There exists a short vector s G M. d such that 

(4) K n 7L d = (s + K) n z d , 

and no facet of any cone s+/Ci, . . . , s + /Q contains any integer points. This gives the disjoint 
irrational decomposition 

JCnz d = (s + /Ci) n z d u • • • u (s + /Q) n z d , 

and so 

i 

(5) a K (x) = Yl x ™ = E^( x ) 

is a rational function. 

For example, suppose that /C is the cone in IR 3 with apex the origin and generators 



wi = (1,0,1), w 2 = (0,1,1), w 3 = (0,-1,1), and w 4 = (-1,0,1) 



(> 
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If we let /Ci be the simple cone with generators wi, W2, W3 and /C2 be the simple cone with 
generators w 2 ,w 3 ,w 4 , then /Q and /C2 decompose /C into simple cones. If s = (|,0, — |), 
then (jlj) holds, and no facet of s + /Q or of s + /C 2 contains any integer points. We display 
these cones, together with their integer points having ^-coordinate 0, 1, or 2. 



s + /Ci 




The cone s + /Ci contains the 5 magenta points shown with positive first coordinate, while 
s + /C 2 contains the other displayed points. Their integer generating functions are 



&s+JC 2 («) 
07c(z) 



x + xz 



and 



Then a s+K;i (x) + cr s+/C2 (z) 



(1 -y*)(l -y l z)(l - xz) 

l + z_ 

(1 - yz){l -y-^z)(l-x-H) 

(l + x)(l-z 2 ) 

(1 — yz)(l — y _1 z)(l — xz)(l — x~ x z) 

a K (x), as 



\x 



xz)(l 



x 



L z) + {l + z)(l-xz) 



x 



xz 



x){\ 



While the cones that appear in the Lawrence- Varchenko formula are all simple, and those 
in Brion's formula are strictly convex, we use yet more general cones in their proof. A 
rational (closed) halfspace is the convex subset of M d defined by 

{x E R d I w • x > b} , 

where w 6 Z d and 6 e 1. Its boundary is the rational hyperplane {x G M. d \ w • x = b}. 
A (closed) cone JC is the interection of finitely many closed halfspaces whose boundary 
hyperp lanes have some point in common. We assume this intersection is irredundant. The 
apex of K, is the intersection of these boundary hyperplanes, which is an affine subspace. 
The generating function for the integer points in K, is the formal Laurent series 



(6) 



K 



E 

m£/C 



X 



This formal series makes sense as a rational function only if K, is strictly convex, that is, if 
its apex is a single point. Otherwise, the apex is a rational affine subspace L, and the cone 
/C is stable under translation by any integer vector w that is parallel to L. If m e /C PI Z d , 
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then the series S%: contains the series 

as a subsum. As this converges only for x — 0, the series Sic converges only for x = 0. 

We relate these formal Laurent series to rational functions. The product of a formal series 
and a polynomial is another formal series. Thus the additive group C[[xf 1 , . . . , x^ 1 ]] of 
formal Laurent series is a module over the ring C[xf x , . . . , x^ 1 } of Laurent polynomials. The 
space PL of polyhedral Laurent series is the Cfx^ 1 , . . . , ccj 1 ] -sub-module of C[[a;f . . . , x^ 1 ]} 
generated by the set of formal series 

{Sic I /C is a simple rational cone} . 

Since any rational cone may be triangulated by simple cones, PL contains the integer gen- 
erating series of all rational cones. 

Let C(xi, . . . ,Xd) be the field of rational functions on C d , which is the quotient field of 
C[xf 1 , . . . , x^ 1 ]. According to Ishida [11], the proof of the following theorem is due to Brion. 

Theorem 7. There is a unique homomorphism of C^f 1 , . . . , x^ l ]-modules 

if : PL — > C(xi, . . . ,x d ) , 

such that <p(Sic) = ctac f or every simple cone JC in M. d . 

Proof. Given a simple rational cone fC — v + (w 1; . . . , w^) with fundamental parallelepiped 
V, we have 

d 

Y[(l-x^)-S K = a v {x). 

i=l 

Hence, for each S E PL, there is a nonzero Laurent polynomial g E C^f 1 , . . . , x^ 1 ] such 
that gS = f E C[xf 1 , . . . , x^ 1 }. If we define f(S) := f/g E C(x\, . . . ,Xd), then f(S) is 
independent of the choice of g. This defines the required homomorphism. □ 

The map (f takes care of the nonconvergence of the generating series Sic when JC is not 
strictly convex. 

Lemma 8. // a rational polyhedral cone tC is not strictly convex, then f(S/c) = 0. 

Proof. Let K, be a rational polyhedral cone that is not strictly convex. Then there is a 
nonzero vector w E 7h d such that w + /C = /C, and so x" ■ S/c = Sic- Thus x™(p(Sic) = (p(Sic)- 
Since 1 — x w is not a zero-divisor in C(xi, . . . , Xd), we conclude that <p(S)c) =0. □ 

We now establish Brion's Formula, first for a simplex, and then use irrational decompo- 
sition for the general case. (A <i-dimensional simplex is the intersection of d+1 halfspaces, 
one for each facet.) 

For a face F of the simplex V, let JCf be the tangent cone to F, which is the intersection 
of the halfspaces corresponding to the d — dim(F) facets containing F . Let be the empty 
face of V, which has dimension —1. Its tangent cone is V. 
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Theorem 9. If P is a simplex, then 

(io) o = ]T(-i) dim ^ F , 

F 

the sum over all faces of P. 

Proof. Consider the coefficient of x m for some m e Z d in the sum on the right. Then m lies 
in the tangent cone JCp to a unique face F of minimal dimension, as P is a simplex. The 
coefficient of x m in the sum becomes 

^(_l)dim(G) _ 

GDF 

But this vanishes, as every interval in the face poset of P is a Boolean lattice. □ 

Now we apply the evaluation map (p of Theorem [7] to the formula (flOl) . Lemma [8] implies 
that ip(S)c F ) = except when F = or F is a vertex, and then tp(S/c F ) = o-/c F (x). This gives 

= -a v (x) + ^ a ic v (x) , 

v a vertex of V 

which is Brion's Formula for simplices. 

Just as for rational cones, every polytope V may be decomposed into simplices Vi, . . . , V\ 
having pairwise disjoint interiors, using only the vertices of V. 

V = Vx U • • • U Vi . 

Then there exists a small real number e > and a short vector s such that if we set 

V 1 := s + (l + e)P and V[ := s + (1 + e)Vi for i = 1, . . . , I , 

then V' n Z d = V fl Z d , and no hyperplane supporting any facet of any simplex V[ meets 7L d . 
If we write /C(Q) W for the tangent cone to a polytope Q at a vertex w, then for v a vertex of 
V with v' = (1 + e)v + s the coresponding vertex of V', we have 1C(V') V > fl Z d = IC(V) V fl Z d 
and so this is an irrational decomposition. Then 

°X(P) v 0r) = a K(V)v{ x ) 

v a vertex of V v a vertex of V 

I 

i=l v a vertex of 7-^ 
I 

= Y av ^ x "> = a -P'( x ) = a r(x) . 

8=1 

The second equality holds because the vertex cones /C("P-) V form an irrational decomposition 
of the vertex cone JC(V') V , and because the same is true for the polytopes. This completes 
our proof of Brion's Formula. 

Consider the quadrilateral Q, which may be triangulated by adding an edge between the 
vertices (2, 0) and (0, 2). Let e = | and s = (— |, — |). Then (1 + e)Q + s has vertices 

(-|,-|), (2,-1), (-1,2 + 1), (4 + 1,2 + 1). 



We display the resulting irrational decomposition. 
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We use the map ip to deduce a very general form of the Lawrence- Varchenko formula. 
Let V be a simple polytope, and for each vertex v of V choose a vector £ v that is not 
perpendicular to any edge direction at v. Form the cone )C^ v V as before. Then we have 

(11) a v (x) = (- 1 ) l ^ (ev)l ^v,v(^)- 

v a vertex of V 

Brion's formula is the special case when each vector £ v points into the interior of the polytope. 
We establish (TTTT) by showing that the sum on the right does not change when any of the 
vectors £ v are rotated. 

Pick a vertex v and vectors £, £' that are not perpendicular to any edge direction at v 
such that £ • w and £ • w' have the same sign for all except one edge direction m at v. Then 
/Q jV and /Q', v are disjoint and their union is the (possibly) half-open cone /C generated by 
the edge directions w at v such that £ • w and £' • w have the same sign, but with apex the 
affine line v + Mm. Thus we have the identity of rational formal series 

S/c 5iV - Sk = -<Sxv, v • 
Applying the evaluation map <p gives 

o-^vO) = -°>c 5 ,, v 0) , 
which proves the claim, and the generalized Lawrence- Varchenko formula (Hip . 



Valuations 

Valuations provide a conceptual approach to these ideas. Once the theory is set up, both 
Brion's Formula and the Lawrence- Varchenko Formula are easy corollaries of duality being 
a valuation. We are indebted to Sasha Barvinok who pointed out this correspondence to the 
second author during a coffee break at the 2005 Park City Mathematical Institute. Let us 
explain. 

Consider the vector space of all functions M. d — > R. Let V be the subspace that is generated 
by indicator functions of polyhedra: 



[V] : x 



1 if x 6 V, 
if x g V. 



We add these functions point- wise. For example, if d = 1, and V = [0,2], Q = [1,3], then 
[P] + [Q] takes the value 1 along [0,1) and (2,3], the value 2 along [1,2], and vanishes 
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everywhere else. 







+ 






' 


L 


2 3 







Already this simple example shows that our generators do not form a basis: they are linearly 
dependent. For V' = [0, 3] and Q! = [1,2], we get the same sum. 



+ 



0123 0123 0123 

But this is the only thing that can happen. 

Theorem 12 (PHI US])- The linear space of relations among the indicator functions [V] of 
convex polyhedra is generated by the relations [V] + [Q] = [V U Q] + [P H Q] where V and Q 
run over polyhedra for which V U Q is convex. 

A valuation is a linear map v: V — > V, where V is some vector space. Some standard 
examples are 



V 


u(P) 


m. d 


vol(P) 


PL 


S v (x) 


C(xi, . . . ,x d ) 


a v (x) 


R d 


1 



That cr-p(x) is a valuation is a deep result of Khovanskii-Pukhlikov [T2] and of Lawrence |14j . 
The last example is called the Euler characteristic. This valuation is surprisingly useful. For 
example, it can be used to prove Theorem [T3l below. 

The most interesting valuation for us comes from the polar construction. The polar V y of 
a polyhedron V is the polyhedron given by 

V y := {x | (x, y) < 1 for all y G V} . 

It is instructive to work through some examples. 



(1) 




The polar of the square 



is the diamond. 



+ 



+ /\ = 
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(2) 





;i,o) 



(0,-1) 

{x | (x, y) < for all y G JC} . 



The polar of a cone K, ... is the cone /C v : = 
(3) Suppose that V is a polytope whose interior contains the orign and ^ is a face of "P. 



Then 




the polar of the tangent cone Kjr 



is the convex hull of the origin 
together with the dual face 
T y := {xeV^ \(x,y) = l}, 
which is a pyramid over T y . 



For this last remark, note that if x G T y and y G JCjr, then (x,y) < (T y ,T) = 1. 
Conversely, if x G /0£-, then (x, .) is maximized over JCjr at T by example (2), and it 
is at most 1 there. 

In these examples, the polar of the polar is the original polyhedron. This happens if and 
only if the original polyhedron contains the origin. 

(4) The polar of the interval [1, 2] is the interval [0, 1/2], but the polar of [0, 1/2] is [0, 2]. 

Now, we come to the main theorem of this section. 

Theorem 13 (Lawrence [2]). The assignment [P] h- > \P W ] defines a valuation. 



This innocent-looking result has powerful consequences. Suppose that V is a polytope 
whose interior contains the orign. Then we can cover V y by pyramids conv(0,^ rV ) over 
the co dimension-one faces T y of V y . The indicator functions of V and the cover differ by 
indicator functions of pyramids of smaller dimension. 



(14) 



[P v ] = y^[conv(0, JF V )] ± lower dimensional pyramids. 



The Euler-Poincare formula for general polytopes organizes this inclusion-exclusion, giving 
the exact expression 



codim.F v +l 



conv(0,J zrV )] . 
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We illustrate this when V is the square. 




If we apply polarity to (1141) , we get the Brianchon-Gram Theorem [HI E] • 

(15) [V] = [/C„] ± tangent cones of faces of positive dimension. 

v vertex 

This is essentially the indicator function version of Theorem [9j but for general polytopes. 
If we now apply the valuation a, and recall that a evaluates to zero on cones that are not 
strictly convex, we obtain Brion's Formula. 

Next, suppose that we are given a generic direction vector £. On a face JF of V, the dot 
product with £ achieves its maximum at a vertex v^J 7 ). For a vertex v of V, we set 

fV(y) ■ = (J relink. 

(The relative interior, relint('P), of a polyhedron V is the topological interior when considered 
as a subspace of its affine hull.) In words, we attach the relative interior of a low-dimensional 
pyramid conv(0, JF V ) to the full- dimensional pyramid conv(0, v v ) which we see when we look 
in the ^-direction from conv(0, JF V ). In this way, we obtain an honest decomposition 

(16) [P v ] = J2 [conv(0,^>))]. 



For the polar of the square, this is 
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To compute the polar of the half-open polyhedron conv(0, J^{v)), we have to write its 
indicator function [conv(0, J~~¥(v))] as a linear combination of indicator functions of (closed) 
polyhedra. If V is a simple polytope, then all the dual faces T y are simplices. It turns out 
that the polar of conv(0, T^{v )) is precisely the forward tangent cone fC^ v at the vertex v. 
So the Lawrence-Varchenko formula is just the polar of (fT6f) . 

This gives a fairly general principle to construct Brion-type formulas: Choose a decompo- 
sition of (the indicator function of) "P v , and then polarize. We invite the reader to set up 
their own equations this way. 

An Application 

Brion's Formula shows that certain data of a polytope — the list of its integer points encoded 
in a generating function — can be reduced to cones. We have already seen how to construct the 
generating function o~ic(x) for a simple cone /C. General cones can be composed from simple 
ones via triangulation and either irrational decomposition or inclusion-exclusion. Given a 
rational polytope V, Brion's Formula allows us to write the possibly huge polynomial o~-p(x) 
as a sum of rational functions, which stem from (triangulations of) the vertex cones. A 
priori it is not clear that this rational-function representation of o~-p(x) is any shorter than 
the original polynomial. That this is indeed possible is due to the signed decomposition 
theorem of Barvinok [I]. 

To state Barvinok's Theorem, we call a rational d-cone K = v + Y2i=i ^>o w « unimodular 
if wi, . . . , W(j G Z d generate the integer lattice Z d . The significance of a unimodular cone K, 
for us is that its fundamental (half-open) parallelepiped contains precisely one integer point 
p, and so the generating function of /C has a very simple and short form 

x p 

0~K.(x) = r ; r . 

v ; (l-x^)---(l-x^) 
In fact, the description length of this is proportional to the description of the cone K,. 

Theorem 17 (Barvinok). For fixed dimension d, the generating function er^ for any rational 
cone K, in IR d can be decomposed into generating functions of unimodular cones in polynomial 
time; that is, there is a polynomial-time algorithm and (polynomially many) unimodular cones 
JCj such that o-)c(x) = J2j e j a iCj ( x )> where ej G {±1}- 
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Here polynomial time refers to the input data of /C, that is, the algorithm runs in time 
polynomial in the input length of, say, the halfspace description of JC. 

Brion's Formula implies that an identical complexity statement can be made about the 
generating function o~-p(x) for any rational polytope V. From here it is a short step (which 
nevertheless needs some justification) to see that one can count integer points in a rational 
polytope in polynomial time. 

We illustrate Barvinok's short signed decomposition for the cone /C := (0, 0) +R> (1, 0) + 
'1,4), ignoring cones of smaller dimension. 
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While K, is the difference of two unimodular cones, it has a unique decomposition as a sum 
of four unimodular cones. 




In general the cone (0, 0) + R>o(l, 0) + R>o(l, n) is the difference of two unimodular cones, 
but it has a unique decomposition into n unimodular cones. 

Arguably the most famous consequence of Barvinok's Theorem applies to Ehrhart quasi- 
polynomials — the counting functions L-p(t) := # (tV H Z d ) in the positve-integer variable t 
for a given rational polytope [I] V. One can show that the generating function J2t>i Lvif) xt 
is a rational function, and Barvinok's Theorem implies that this rational function can be 
computed in polynomial time. Barvinok's algorithm has been implemented in the software 
packages barvinok [T7] and LattE [8|. The method of irrational decomposition has also been 
implemented in LattE, considerably improving its performance |13j . 
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